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We study the properties of the Fraunhofer diffraction patterns produced by Gaussian beams 
crossing spiral phase plates. We show, both analytically and numerically, that off-axis displacements 
of the input beam produce asymmetric diffraction patterns. The intensity profile along the direction 
■ of maximum asymmetry shows two different peaks. We find that the intensity ratio between these 

' two peaks decreases exponentially with the off-axis displacement of the incident beam, the decay 

, being steeper for higher strengths of the optical singularity of the spiral phase plate. We analyze 

^N) ■ how this intensity ratio can be used to measure small misalignments of the input beam with a very 

^ ' high precision. 
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I. INTRODUCTION 



Oh] Optical vortices (OVs) appear in light beams carrying screw wavefront dislocations (vortex beams) [l|. The surface 
O ,: of constant phase of a vortex beam has an helical structure and presents phase singularities endowed with topological 
c/3 ■ charge. Beams harboring OVs carry also a quantity of orbital angular momentum (0AM) [2| associated to the 
^ , precession of the Poynting vector around the vortex axis 0] . 

^ ' OVs have attracted an increasing interest in applied physics [1, d, [1, 0) B and also for astronomical applications 
[3 EH, [H E B [3 • In fact , they can be easily produced in light beams with the help of specific optical devices that 
have a central optical singularity. Among these optical elements, the most efficient ones are fork holograms (FHs) 
and spiral phase plates (SPPs) [l3| ■ Laguerre-Gaussian (L-G) modes have been often used to describe the beams 
produced with such devices. However, a more precise description of the diffraction patterns produced by an SPP [l8| 
or a FH [19, 20] is provided by hypergeometric (Kummer) functions. We shall use this approach in this Paper. 

Consider an input beam with an amplitude distribution symmetric about the propagation axis. When such a beam 
intersects an SPP or a FH perpendicularly and exactly on-axis, it produces a circularly symmetric beam with a central 
: dark region, where the field amplitude is zero. Any misalignment with respect to the central discontinuity would then 
■r^lj- ■ produce an asymmetry of the observed intensity distribution [21,] and the topological charge of the correspondent 
[ off-axis OV may have a non-integer value This changes also the 0AM originally carried by the beam [23|, thus 
• producing an 0AM spectrum [2J, |2^ . 
^\ ' The sensitivity of a vortex beam to displacements of the input beam has been proposed as an indicator of nanometric 
. shifts in a speckle pattern [26j or to be used as a non-interferometric method for the correction of small surface 
^ ' deviations on spatial light modulators 12^ . Similar results could be obtained from the analysis of the mean square 
. ^ , value of the resulting 0AM spectrum [28| . We also proposed a method to measure very small displacements based 
on the degree of asymmetry of the intensity pattern of an off-axis vortex beam ^13, .29]. However, detailed analytical 
H ] studies of the actual structure of off-axis OV produced with SPPs of FHs were initiated only very recently [30] . In 
■ - - 1 this Paper, we extend the analysis of the intensity distribution of an off-axis vortex beam generated under Fraunhofer 
diffraction conditions. We then derive a more convenient formalism of our method for the detection of off-axis 
displacements, in view of future applications with optical imaging devices. 
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FIG. 1: The geometrical configuration adopted to study the Fraunhofer diffraction of a Gaussian beam beyond a spiral phase 
plate (see text for details). 



II. FRAUNHOFER DIFFRACTION OF A GAUSSIAN BEAM INTERSECTING A SPIRAL PHASE 

PLATE ON-AXIS 

In this Section we revisit the problem of the Fraunhofer diffraction of a monochromatic Gaussian beam intersecting 
an SPP exactly on-axis [i^. The geometrical configuration adopted here is sketched in Fig. [T] The SPP is placed in 
the {x,y) plane and the central optical singularity coincides with the origin of the {x,y) coordinate system. To take 
advantage of the circular symmetry of the geometry, in the following wc will use the circular coordinates (r, 6) defined 
by X ^ rcos6 and y = rsm6. Thus, the transmission function of the SPP is a complex function of the azimuthal 
angle 9: 

Teie)=e''\ (I) 

where £ represents the strength of the optical singularity. We consider a Gaussian beam propagating along the z axis 
(that is also its symmetry axis) perpendicular to the SPP plane. This choice is important for practical applications, 
i.e. laser beams or starlight beams affected by atmospheric turbulence. We then assume that the field amplitude 
distribution in the SPP plane is 

Acir) = ce-'^'/'"' , (2) 

where c is a complex factor, eventually dependent on the z coordinate, and w is related to the beam width. The 
observation plane {x', y') is located at a distance / ^ oo beyond the SPP or, better, at the focal plane of a lens placed 
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just behind the SPP (in this case, / would coincide with the focal length of the lens). The scalar field of the beam in 
the observation plane can be obtained from the Fourier transform of the product between the amplitude of the input 
beam and the SPP transmission function: 

(3) 



= const J J AG{r)Te{e) 6^'''''^''°''^^-''^^/^ rdrd0 



where k = 2tt/X is the wave vector, A is the wavelength and (p, /?) are the circular coordinates in the observation 
plane. To simplify the calculations, the SPP is considered as infinitely extended in the {x,y) plane. In addition, we 
use the scaled radial coordinate r' — kr/f and introduce the quantities = f /{kw) and cq, the latter containing all 
the multiplicative constants. In this way, Eq. [3] becomes 



ueip,P) = Co 



2tt 



(4) 



The integral involving the angular coordinate 9 can be evaluated by using the definition of the Bessel function of the 
first kind J„(z). Thus, the integral involving the spatial coordinate becomes a particular case of the Weber-Sonine 
formula [31]. The final result of the integration can be expressed in terms of the confluent hypergeometric function 
of the first kind. However, a more useful expression of the amplitude distribution of the output beam is obtained by 
using the modified Bessel function of the first kind I,y(z). By introducing the quantity rj = p/[2wq), the final result is 

M 



ue{p,P) = coi 



2wl 



(5) 



The presence of the phase factor e'^'^ implies that the output beam has an ^-charged OV nested inside. 

Beams of this kind, also known as 'Kummer beams' (20l |. are different from the commonly used Laguerre-Gaussian 
(L-G) beams [s^, HI]. If £ = 0, the Bessel functions of half-integer index in Eq. [5] can be expressed in terms of the 
hyperbolic functions and combined together to give an exponential. In this case, the amplitude distribution of the 
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output beam is still Gaussian, i.e. uq ~ e^'' . If £ 7^ 0, we can derive an useful approximation for 77 ^ by using the 
series expansion of 1,^(2;) [sil ]: 



M.)^(l)' 



y- (^/2) 
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(6) 



We can recognize that, near the z axis, the amplitude of an on-axis Kummer beam carrying an OV with topological 
charge £ could be represented by a superposition of amplitudes of L-G modes with p = 0: 
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(7) 



The dominant term is represented by an L-G mode with index £, while higher order terms are L-G modes with indices 
e + 2m{m = 1,2,...). 

A. Properties of the intensity distribution 

The intensity distribution of an on-axis Kummer beam is axially symmetric around the z axis and is described by: 



h{p, 13) ~ \ue{p, P)\ =Con 



4wt 



I.-. ( ^ 



(8) 



As for L-G modes, the intensity pattern of a Kummer beam has an annular shape, with a central dark region where 
the intensity is zero. However, there are some fundamental differences between the two analytical descriptions. For 
a Kummer beam the behavior of the intensity at large distances from the z axis is ^ ?7~^, while for an L-G mode 
it decreases exponentially. Moreover, the radius of maximum intensity of an L-G mode is pmax ~ •\/-^/2, where the 
intensity attains the value /(pmax) ~ i^e~^/£l, while for a Kummer beam pmax is found by numerically solving the 
transcendental equation 



(9) 




FIG. 2: Properties of the intensity distribution of on-axis Kummer beams having ^ = 0, 1, . . . , 10. (a) Plot of the radius of 
maximum intensity pmax (in units of wo) vs. £. (b) Plot of the intensity calculated at pmax for different values of £. The 
maximum intensity for 1 = has been set to unity. 



The calculation of the radii of maximum intensity obtained for a set of values of the topological charge ^ = 0, 1, . . . , 10 
(see Fig. [5^) suggests that pmax is linearly dependent on £: 

= (0.37 ± 0.01) + (0.470 ± 0.002)£ . (10) 

A similar result was found also for OVs produced by a plane wave intersecting a finite circular phase mask 6]. Fig. 
shows the intensity calculated at pmax for the same set of topological charges. 



B. The effects of off-axis displacements 



When the input Gaussian beam is displaced off-axis, so that its symmetry axis does not coincide with the z axis of 
Fig. [H the intensity pattern produced in the observation plane is modified. The misalignment of the input beam can 
be decomposed into a translation in the {x, y) plane and an inclination angle to with respect to the z axis. However, 
if Lo is small, the modifications induced in the intensity pattern of the output beam are negligible [30| . For this 
reason, in our calculations we will consider only lateral displacements of the incident beam. Let us then assume that 
the Gaussian beam intersects the SPP in the position (rofi,6'off), as shown in Fig. [1] The field of the output beam 
resembles that of Eq. [5] obtained under on-axis conditions: 
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FIG. 3: Example of the far-field intensity pattern of a numerically-simulated off-axis vortex beam produced by an ^ = 2 SPP. 
The two diamonds indicate the positions of the two different intensity peaks A and B. (a) Contour plot of the intensity 
distribution obtained in the observation plane, (b) Intensity profile of the off-axis beam extracted along the x' direction. A is 
the difference between the intensities calculated in B and A. Examples of error bars are reported for both the intensity peaks 
(see text). 



Here, r^g is the scaled radial coordinate obtained from Toff, while the quantities 7 and t/j are defined as [31 



In this case, the additional exponential factor and the complex value of 7^ produces a phase singularity which is located 
neither on the beam axis, nor in the origin of the (x', y') plane, but shifted in a position (p, (3) — i^w^r'^^, OqH + 7r/2) . 
As a result, the intensity distribution of the output beam becomes asymmetric [34], showing two different peaks along 
the direction of the vortex core in the (a;', y') plane. Fig. [3^ shows an example of an off-axis OV produced with an 
t = 2 SPP. The lower and the higher peaks are labeled with A and respectively. 

Now, since it is difficult to find analytical solutions of Eqs. Illl and fT^ we decided to perform numerical simulations. 
The width w of the input Gaussian beam was parametrized in function of the full width at half of the intensity 
maximum, 2a, such that a = Wi/ln 2/2. We used values of the topological charge induced by the SPP in the range 
£ = 0, 1, . . . , 10, since with higher values we get misleading results using the two-dimensional Fast Fourier Transform 
algorithm. For each £, we considered a number of off-axis displacements rog/a of the input beam ranging from to 1 
and computed the intensity patterns of the resulting beams. 

We checked the consistency of our numerical simulations by comparing them to the analytical models (Eq. Ilip 
for a number of values of i and off-axis positions. To this aim, we previously normalized the intensities of both the 
simulated and the analytical patterns to the corresponding maximum values. Therefore, the B peak always has a 
normalized intensity equal to one (obviously, both the peaks A and B will have the same unity intensity if roff/a ~ 0). 
The residuals of the subtraction of the theoretical intensity patterns from the simulated ones are typically within 
10"'* for positions close to peaks A and B. We will then assume this quantity as the intrinsic error of our numerical 
simulations. 

For all the numerically simulated OVs we obtained the intensity values at the two peaks and calculated the quantity 
R defined as the ratio between the intensity I a of the lower peak and the intensity Ib of the higher peak. We find 
that R rapidly decreases as the off-axis displacement increases for all the topological charges considered. The graphs 
showing the dependence of R on ros/a for £=1,2, 3, 4, 5 are plotted in Fig. H) All the curves are well represented by 




(12) 
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FIG. 4: Plot of the peaks intensity ratio R vs. the off-axis displacement of the input Gaussian beam obtained for different 
values of the topological charge induced by the SPP. Horizontal lines are drawn at the maximum values of R detectable at the 
la confidence level for 3 values of the SNR ratio of the B peak. 



a simple exponential function: 

i? = fcie-'=^''°«/'' , (13) 

where parameters ki and k2, obtained by best fitting the simulated curves, are listed in Table HI From these results, 
it appears that ^2 depends on the topological charge £ as 

fc2 = (4.64 ± 0.05) - (2.9 ± o.2)e-(°-^^±"-°5)^ , (14) 

while ki seems to remain equal to unity. 



III. THE SENSITIVITY OF THE METHOD TO REVEAL OFF-AXIS DISPLACEMENTS 



The off-axis displacement of the input Gaussian beam with respect to the central singularity of an SPP results in 
an asymmetry of the far- field intensity pattern. Eq. [T3] reveals that the parameter R represents an extremely sensitive 
tool to detect such very small displacements. 
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TABLE I: Best-fit values of parameters fci and k2 in Eq. [13] obtained from least-square fits of the simulated curves of Fig. [4] 
The associated errors are given at the la confidence level. 



e ki k2 

T 1.002 ±0.001 2.808 ± 0.006 

2 1.0012 ± 0.0007 3.534 ± 0.004 

3 1.0004 ± 0.0005 3.928 ± 0.003 

4 0.9995 ± 0.0008 4.259 ± 0.005 

5 0.9990 ± 0.0009 4.338 ± 0.006 

6 0.999 ± 0.001 4.432 ± 0.008 

7 0.998 ± 0.001 4.456 ± 0.008 

8 0.998 ± 0.001 4.509 ± 0.009 

9 0.998 ±0.001 4.573 ±0.01 

10 0.998 ± 0.002 4.756 ± 0.02 



Let us suppose to have an input Gaussian beam which symmetry axis is perpendicular to the surface of an SPP and 
observe the correspondent far-field image with a photoelectric detector like a CCD. In this way, if the beam is displaced 
off-axis, we can measure the intensity ratio R of the two different peaks, as defined in Sect. IIIBI The precision of 
this measurement is mainly limited by the signal-to-noise ratio (SNR) achieved in the observation, while additional 
errors might be introduced by construction imperfections of the SPP. The latter issue results in inhomogeneities of the 
observed intensity distribution. For this reason, efforts are currently made to improve the production quality of SPPs 
[S^ll^ls^]- However, here we assume to use an ideal SPP so that the only limitations are due to the SNR. The noise 
of a CCD detector is mainly represented by the photon shot noise [111. Assuming a pure Poissonian distribution of 
the collected photons, the uncertainties associated to the intensities Ia and Ib of the two peaks can be approximated 
by the square root of the signals, i.e. cta = a/Ta and as = \/Ib- We may then recognize Ia 7^ Ib at the na confidence 
level when A = Ib — Ia > noA ± naB (see Fig. [5]d), that means 



A>n[yiA + ^/lB] . (15) 

If we introduce parameter i?, this equation can be rewritten as a function of the only SNR associated to the intensity 
of the highest peak (SNR(i3)). We find that the maximum peaks intensity ratio measurable at the na confidence level 



As useful examples, in Fig. [4] we draw three horizontal lines corresponding to the maximum R values 0.44, 0.81 and 
0.98 detectable at the la level for SNR(B) = 3, 10 and 100, respectively. By combining Eq. [Hlwith Eq. [H we 
finally obtain the expression for the minimum off-axis displacement detectable at the na confidence level: 



^>-^ln 
a k2 



Yiy^ SNR(B) 



(17) 



One general outcome is that, for a fixed SNR(i?), OVs with higher £ values allow the detection of smaller off-axis 
displacements. This effect is more significant at low SNR regimes, when the maximum measurable R is small and the 
curves in Fig. [Jare more spatially separated. As SNR(i3) increases, the advantage obtained by using high values of the 
topological charge becomes negligible. In fact, if we assume SNR(i3) above 10, we might reveal off-axis displacements 
< 0.1a for all t values. Instead, considering the lowest acceptable value SNR(_B) — 3 for signal detection, we can 
detect off-axis displacements of ^ 0.3a for £ — 1 at the la confidence level. 



IV. CONCLUSIONS 



In this Paper we have analyzed the properties of the Fraunhofer diffraction pattern produced by a Gaussian light 
beam crossing an SPP. When the input beam is perfectly aligned with the central singularity of the SPP, the resulting 
beam is a Kummer beam with a symmetric annular intensity distribution. Instead, an off-axis displacement of the 
input beam produces an asymmetry in the far-field intensity pattern. In particular, the intensity profile along the 
direction of maximum asymmetry shows two different peaks. We have found that, for all the values of the topological 
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charge considered, the ratio R of their intensities changes exponentially with the ofF-axis displacement of the input 
beam. We have quantitatively analyzed how the SNR associated to the highest peak affects the sensitivity of the 
ratio R in revealing very small misalignments of the input beam. In particular, we have found that higher values of 
the topological charge £ generally provide better resolutions, especially for low SNR regimes. We suggest that this 
method could find interesting applications in high-precision positioning systems. Note that similar results can be 
obtained also by using other spatial properties of transverse laser modes 39] . Finally, the sensitivity of OVs could be 
used in astrometry, by placing an SPP at the focal plane of a telescope. 
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